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ABSTRACT 


Stability in the sense of Lyapunov is investigated for systems 
of ordinary differential equations. Various established results 
are discussed and a new set of sufficient conditions for a stability 
comparison theorem is obtained. The treatment is principally 
concerned with linear systems, though some results for non-linear 


systems are considered briefly. 
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I. INTRODUCTION 


Stability was originally a pnysical concept, believed first to have 
been considered by Lagrange in connection with the equilibrium states 
of conservative systems. Equilibrium states are those in which the 
velocities and accelerations are zero, describina a state of rest. 

A svstem in equilibrium may react in one of two manners when influ- 
enced externally; it may exhibit motions which remain in a small 
reaion about the equilibrium position (which does not preclude the 
nossibilitv that the system return to its state of equilibrium), 

or it may exhibit motions which cause the system to increasinalv 
deviate from its equilibrium state, never returning to it. In the 
former case, the equilibrium state is said to be stable, while the 
latter case describes an unstable state of equilibrium. 

The above remarks are, of course, not rigorous and are intended 
solely as an intuitive guide to the meaning of stability of equili- 
brium states of nhvsical systems. Physical systems are almost 
always described by a system of differential equations and in 
Studving the behavior of a physical system, we ultimately study the 
behavior of the solutions to the differential equations which 
describe the system. Tnus the study of stability of a physical 
system is reduced to a consideration of the stability of solutions 
of differential equations, and stability theorv is todav an intearal 
part of that segment of analysis known as the theory of ordinary 


differential equations. 


Stability theory is quite advanced todav and there are numerous 
concepts and approaches which vary with both the form and order 
of the system of differential equations. For example, the Poincare- 
Bendixson theory is an analytical-topoloqical apnroach which is 
particularly useful for second order autonomous svstems, while for 
non-autonomous systems, the concept of stabilitv in the sense of 
Lvanunov is utilized and the approach is purely analytical. This 
paper is concerned with a general study of stability in the sense of 
Lyaounov, with particular emphasis on systems of linear differential 
equations though some results for nonlinear systems are considered. 
Section II deals with a review of familiar facts from the theory of 
linear ordinary differential equations. Stabilitv and boundedness for 
linear systems are discussed in section III. Non-homoqeneous svstems 
are considered in section IV, and tne concents of uniform and restric- 
tive stability are treated in section V. Some results for non-linear 
systems are given in section VI, and concluding remarks follow in 


section VII. 


IIT. PRELIMINARY REMARKS 


We consider the behavior of solutions of the nth order linear 
non-homogeneous system 
(2.1) x' = A(t)x 
where x = x(t) is an n by 1 column vector and A(t) is an n by n matrix 
with elements a; ;(t), i,j = 1,...n,. which are real valued functions 
continuous on the half-line [t, 2°). 

The norm of x, A is defined by 
(P27) || x(t) | | 
22211) | A(t) | | 


E(x. (t) | 
EIA, 5 (t)| 
and note that both ||x|]| and |]A|| are functions of t. The usual 


properties of norms are satisfied and are stated here without proof. 


For ||x|| they are, 





CAR (STN) eigen 
(2.341) LHixt] = Hyll{ < Pbeeyll <tletl + Tyll. 
(Dash) Jlax}|{ = Jol | Ix] for a a real number, 
(Za a learece ts aecorilaloc( €) |i) city, 

and for ||Al], 

41.) tee Oemmanae| (Alle = Oni ff A(t) = 0, 
(2.473) Meese sal) + | )Bi 
(2.4473) [JAB] | < TIAT| [{BII; 

(2.4iv) HIAx|] < TAI] IxI1, 

(2.4v) ||oA]} = Jo} ||A]| for o a real number. 


Canital letters will denote n by n matrices, lower case letters 


n by 1 column vectors, and Greek letters real numbers unless otherwise 


noted. Primes (') denote differentiation with resvect to t; al denotes 
the transpose of A; tr(A) denotes the trace of A; det A denotes the 
determinant of A; and adj A decotes the adjoint of A. Lower subscripts 
denote components of a vector or matrix, for example, Xs is the ith 
component of the vector x. 

The questions of existence and uniqueness of solutions to the 
system (2.1) may be dealt with quite briefly. Solutions to (2.1) 
satisfying the initial condition x(t.) exe exist and are unique on 
any interval in which A(t) is continuous. Hereafter, continuity of 
A(t) on [t, 2°) will be assumed so that existence and uniqueness of 
solutions for systems of this form are always assured. 

Linear homogeneous systems of differential equations are charac- 
terized by the fact tnat any linear combination of solutions is also 
a solution. The space of all possible solutions to the system (2.1) 
constitutes an n-dimensional vector space. Thus there exists n 
linearly independent solutions, nl ee we which form a basis for 
this vector space of solutions. Any particular solution x(t) of (2.1) 
can then be expressed as a linear combination of elements in this 
basis. Since any n linearly independent solutions constitute a basis, 
it is convenient to pick a particular set and define the fundamental 
matrix X(t) as follows. 

(2.5) K(t) = o4(t) (i,j = Ty..eon), 

i.e., the columns of X(t) are n linearly indenendent solutions of 
(2.1) and these n solutions 4! (i = 1,...,n) are chosen to satisfy the 
initial condition 63(ty) See eee, ih) en eal X(t.) = ee 


1J 
the identity matrix. 


From (2.5), we note that X(t) satisfies the matrix differential 
equation 
(2.6) ee EK, 
and that any solution x(t3t, sx) of 297) Satistying™the™ initial 
condition x(t) = xX) may be written 
(2.7) wale) = X(t)x,. 
From the formula for the derivative of the determinant of a 
matrix, 


d det x(t) = tr [A(t)] det x(t), 
t 


d 
we obtain the Jacobi-Liouville formula 
1E 
(2.8) det X(t) = det X(t.) exp fs tr A(s) ds 
ty 
t 
= exon fs tr A(s) ds, 
iB 
O 


and note from (2.8) that X(t) is non-singular for t a 
Any non-singular matrix W(t) which is a solution of (2.6) mav 
serve as a basis matrix though if W(t) # I, weno longer have (2.7) 


but instead 


where ae # Xo 
From tne formula for the derivative of the inverse of a matrix, 
a. Ux! (t)] = =x '(t) x(t) xT'(E), 
we obtain the result that if Y(t) = x Neeyl then Y(t) satisfies 
the matrix differential equation 


T 


(2.9) ie =e -AGe)” Ye 


and since Y(t) is non-singular and Y(t.) = I, we note that Y(t) 
is the fundamental matrix to the adjoint system of (2.1) qiven by 
(2.10) xt = A(t)! x, 

The fundamental matrix X(t) provides an explicit form for the 
solution of a non-homogeneous system of the form 
(Zale) x' = At) xeeesbttt) 
as the following proposition illustrates. 


Proposition 2.1. The solution x(tst. 4x of the system (2.11), 


~ 
where A(t) and b(t) are continuous on Lta se) may be written 
8 ig eee eee) a! ee) ee bisaias, 
t 
0 
where X(t) is the fundamental matrix of the homogeneous system 
Xe ee oe 
Proof: Direct substitution. 
A similar formula may be obtained for non-linear mixed systems 
of the form 
(are) CARA Ke, AEE E 
Proposition 2.2. Let f{t,x) be continuous for t cil and 
||x|| < €, and A(t) continuous for t > t,. Then the solution 
x(t; t_,x_) of (2.12) may be written 
O20 t 


x(t) = X(t) x +f X(t) X-"(s) fls,x(s)) ds, 
ic 


0 
where X(t) is the fundamental matrix of x' = A(t) x. 
Proof: Direct substitution. 
Both Propositions 2.1 and 2.2 simplify in the case A(t) = A, 


a constant matrix and ty = 0. For then, the matrix X(t) x7 05) 


10 


satisfies (2.6) with initial condition 
Me) Xs) = 
at t = s, and the same follows for X(t-s). Hence by uniqueness, 


(2.13) NGpieeas) = (tes), 


for A a constant matrix and X(0) = I. 
Other well known results from the theory of ordinary differential 


equations will be stated when needed. 


| 


TTI. LINEAR HOMOGENEOUS SYSTEMS 

Consider a general first order system of ordinary differential 
equations of the form 
(Sia) ole earn tea) 
where x(t) is ann by 1 column vector and f(t,x) is ann by 1 
column vector defined and continuous in the region I’ = {(t,x): 
O<t<m, ||x]| < ot. 

We wish to define a concept of stability for this system. One 
definition which immediately comes to mind miaht be to define the 
system to be stable if it exhibits only finite motions, i.e., all 
solutions remain bounded as t-~-. However, sucn a concept is not 
particularly useful as it distinguishes variational effects only as 
to their being finite of not. For example, an earthquake would be 
Stable in this sense. 

What we wish to embody in a stability concent is a means bv which 
we can measure the effects of small variations on the system (or a 
particular solution of the system). One criterion we would desire of 
a stable system is that small variations not affect the system in a 
drastic manner. This is the orimary consideration embodied in the 
concept of Lyapunov stability, which imposes stringent restrictions 
on the solutions of (3.1), by requiring that solutions which are 
initially close together remain close together for all future time. 


Formally, we define stability as follows. 


12 


Definition 3.1. Let x(tst x) be a solution of (3.1) satisfying 


(3.21) Gey exists for © > tos 
and 
eecin) (Gaither 6©6 for t > t- 


Then x(t) is said to be stable (in the sense of Lyanunov) if 
(Sect) there exists a>0 such that every solution 
y(t3t sy.) Exists on [t»~) and (t,y(t)) 
whenever ||x. - y || < as 
and 
(3.2iv) given e>0 there exists 6 = é(e,f,x), 
O<d<a, such that ||x(t) - y(t)|| <«€ for t ae 
for all solutions y(t) with 11x, - YI | < 6. 

Strictly speaking, the above defines stability at the right, since 
the behavior as te» is considered; replacing the interval [t»~) by 
the interval (ocr would define stability at the left. Note 
from (3.2iv) that stability is equivalent to uniform continuity of 
x(t3t..x,) on [t,»°) with respect to the initial vector Xo" A solution 
which is not stable is said to be unstable. Hereafter, stability will 
always refer to Lyapunov stability at the right as defined above. 

A special case of stability occurs when solutions which are 
initially close approach one another as t>~. This leads to the 
following definition. 

Definition 3.2. A solution x(tst..x,) of (3.1) is said to be 
asymptotically stable on [t,»°) if it is stable and in addition there 
exists 6' > 0 such that |.|x(t) - y(t)||+0 as t+~ whenever yv(t;t.y,) 


° * ° ’ 
is a solution with IIx, - spall ausite. 


(1s: 


Note that by definition, a solution which is asymptotically stable 
is stable. If all solutions of the system (3.1) are (asymntotically) 
stable, then we say that the system is (asymptoticallv) stable, while 
if there exists an unstable solution to (3.1), we will call the system 
unstable. 

The following examples illustrate some of these concents. 

Examnole 3.1. Consider the one-dimensional system 

x' = =x 
on [0,0). The general solution x(t30,x)) is given by x(t) = Xy @ is 
Som una 
=e 


||x(t) - y(t)|] = |x, - yol e 


> 40 Sy 


| — 
Conditions (3.2 i,ii, and iii) are satisfied as is (3.2iv) with 
6 = ¢. Thus the system is stable. 
In addition, 
i}x(t) - y(t)]] = [x e: > diiata cies 


so that the system is asymptotically stable. 


Oo. vi 


Example 3.2. Every non-trivial solution of the one-dimensional 
System 
eee 0) royal. [Oren 
is of the form x(t30,x) = and thus 
HIx(t) - y(t)|| = |x, + y 


Conditions. (3.21,i171, and iii) are satisfied as isml@e2iy eo aoineuatD 


o| - 


Oe= 6. 
However, 


lim ||x(t) - y(t)|] = |x, -y | > 0, 


to 


so that no solution is asymptotically stable. Thus the system is 


stable. 
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Example 3.3. Non-trivial solutions of the one-dimensional system 


a= Ome Onc) 
are of the form x(t30,x,) = xe" and thus 
_ c 
[Ix(t) - y(t)|] = |x, - yp le 


which can be made arbitrarily large independent of Xo - Yo! by taking 
t large enough. Thus the system is unstable. 
Examnle 3.4. Cesari [2] shows that the trivial solution z(t) = 0 
of the one dimensional system 
x' = ie on [0,-) 
is not stable since every (real) solution x(t30,x,) is of the form 


x(t) = 0 


1-tx, 


and solutions with Xen 2 0 cease to exist at t = (x,)7 in violation 
Oi ad sae igi) 
For linear systems of the form 
ioe 3) x = Bit) x + bigs 
with A(t) and b(t) continuous on [t 22); (3.2i1,i17, and iii) are always 
satisfied so that (3.2iv) becomes the determining factor as to whether 
or not solutions to these systems are stable. For these systems, 
we have the following useful proposition. 
Proposition 3.1. If the linear system (3.3) has a stable solution 
then all its solutions are stable, i.e., the svstem is stable. 
Proof: Let x(t3t sx,) be a stable solution and y(tst sy) 
be any other solution. Let « > Q. 
Ie z(t3t,»Z,) is a solution of (3.3), then x(t) - y(t) + z(t) 
is also a solution since the system is linear. Then since x(t) 


is stable, there exists 6 > 0 such that if 


[ts 


ibe) Loot acamanlatees acta) ie =i | ly onal 


mr 
then 
}[x(t) - Ex(t) - y(t) + z(t)J|| = [ly(t) - z(t)]| 
< € 

which says that y(t) is stable. Thus all solutions are stable. 

Corollary |... ‘Propositron 37) rémaias™truemif stable ism@meplaced 
by asymptotically stable. 

Proof: The proof is identical to that above. 

We next define the concept of boundedness of solutions to 3.1. 

Definition 3.3. A solution x(t) of (3.1) is said to be bounded 
if there exists M > 0 such that || x(t)|| <M, t > to: 
i.e., if ||x(t)]] is a bounded function of t on [t,.~). aH | 
solutions to (3.1) are bounded, we say that the system is bounded. 

Stability and boundedness are independent concepts as the 
following examples illustrate. 

Example 3.5. All solutions of the one dimensional system 

me = 92 on | One) 

are of the form x(t3t x) = es t¢ and are stable yet unbounded 
One Olen) 

Example 3.6. Cesari [2] shows that every non-trivial solution of 


the two dimensional system 


ay 

y' = - 1/2 x (x¢ + a. + i ae 
or equivalently, 

2) 2/2 aaa 
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is of the form 
aay a= Cc Sin (cuta@e ‘c # O 

and is thus bounded; none of these solutions is stable. 

For linear homogeneous systems of the form 
(G4) x = Wee) x 
stability and boundedness are equivalent as the following theorem 
illustrates. 

Theorem 3.1. The system (3.4) is stable iff it is bounded. 

Proof: Assume the system is stable. To show that it is bounded, 
it suffices to show that there exists M > 0 such that ||X(t)|| < M 
where X(t) is the fundamental matrix of the system. Since the system 
is stable, in particular, the trivial solution z(t) = 0 is stable. 
Thus if « > 0 there exists § > 0 such that every solution x(t3t,,x,) 
with ||x || < 6 satisfies ||x(t)|| < « for t>t,. In particular, let 
x(t) be that solution with initial vector xe given by 

6 1 


TG 


where e| is the vector whose jth component is given by e, = Si: 


Then ||x,|| < 6 so that 


Pee (Dix |Ih<uc, t >t. 
But for this choice of Xo? X(t) Xn > yl where al is the ith column 


of X(t). Thus 


slit teat . 
a = 0 


This procedure is valid for 1=1,...,n so that 





L1x(t) || Hol(t)|| < Fe = m 


i 
a) a) 
af 


l= 


and hence X(t) is bounded. Thus the system is bounded 


17 


Conversely, assume the system is bounded so that X(t) < M, 
fOReE > t- Then if X(t3t x) and y(tst,sy,) are two solutions to 


(3.4) 


||x(t) - y(t)] | AG ee oy) |i] SM) 


so that if ||x, - yoIl < mw > then ||x(t) - y(t)|| <e fort aie 


si 


which shows that the system is stable. 

Corollary 1. The system x' = A(t) x is asymptotically stable 
iff X(t)-0 as tem, where X(t) is the fundamental matrix of the system. 

Proof: The reverse implication is obvious. To show the forward 
implication, note that the hypotheses imply that the trivial solution 
is asymptotically stable and the result then follows as in the proof 
of the first half of Theorem 3.1. 

Before proceeding to consider various results concerning stability 
in linear homogeneous systems, it is convenient to introduce the 
following lemma, known as Gronwall's Lemma, which is an extremely 
useful tool in stability theory. 

Lemma 3.1. If f(t) and g(t) are real valued non-negative functions, 
f(t) is continuous on [t»@) and g(t) is integrable in everv finite 


interval contained in pts?) and if for some constant a > 0 


c 
(3.5) Tl Cymer 1S gic) cram 
Ee 
then 
t 
(3.6) f(t) < a exp f Giciieasr 
t 
0 


Proof: If a> 0, then we have 


f(t) g(t) sgt). 


18 


Integrating the above between the limits to and t yields 


1 ie 
In [a + s f(s)g(s)ds] - Ina < fs g(s)ds, 


LG 0 
or taking exponents, 
iG 1p 
ar fomeoygrsjas <a exp fagis\de, 
et c 
0 0 


and thus 


u 
f(t) <a exp fs g(s)ds. 
tO 
If ao = O, then (3.5) holds for all real numbers g > 0, and thus 


(3.6) holds which implies as g+0 that f(t) = 0. This completes the 
proof. 

Utilizing Lemma 3.1, it is possible to prove the following 
proposition. 


Proposition 3.2. If sj{|{A(t)|| dt < , then the system x' = A(t) x 
t, 
Ws" Cable. 


Proof: Let X(t) be the fundamental matrix of the system, Then 


from (2.6) 
ie lr ff Mic) X(S) OS, 
so that 


fe 
J[XCt)}f sn + £  [fACs)|] [[Xs)|] ds. 
by 
Applying Gronwall's Lemma yields 


i 00 
}|X(t){]| < nexn s |fA(s)|{| ds < nexp s ||A(s)|| ds 
A t, 
and thus X(t) is bounded. By Theorem 3.1, the system is stable. 
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We next consider the linear homogeneous system with constant 
coefficients 
(3.7) x =e 
whose fundamental matrix X(t) is given by 
(3.8) X(t) = exp [(t-t,)Al, 


where exn (A) is defined by 


exp (A) = I+. ¢ a ne 


n=] 

The stability of the system (3.7) is directly related to the nature 
of the roots of the characteristic polynomial of A given by 
(3.9) det (A - x1) = O. 
Since (3.9) is an nth degree polynomial equation in , there exists 
n roots to this equation (including repeated roots) which may be 
real or complex. Corresponding to each root is a solution of (3.7) 
and the n roots thus furnish n linearly independent solutions to 
(3.7). If A is a real root of multiplicity k, then corresponding 
to » are the k linearly independent solutions 

ert. ter. — pk- Tat 
while if Atip is a complex root of multiplicity m, then so is 
A-ip and corresponding to these 2m roots are the 2m linearly indeven- 
dent solutions 
et eos ut, etc in ut; ter *cos oe te ’sin ire ae Ml hos Wh, el sin ut. 
Since sin t and cos t are bounded, and since for a, 8 >0, 
yt 


lim te = 0, we note that if the real parts of the eigenvalues of 


t->0o 
A (roots to (3.9)) are negative, then all solutions of (3.7) will 


converge to zero. When this is the case, we have the following 


definition, 
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Definition 3.4. If the real parts of the eigenvalues of A are 
negative, then we say that the characteristic polynomial of A is 
stable. 

When the characteristic polynomial of A is stable, then the 
solutions of (3.7) are of exponential order as the following lemma 
illustrates. 

Lemma 3.2. If the characteristic nolynomial of A if stable, 


then all solutions of (3.7) satisfy 


(3.10) Wectip DN sie - at acs rey, 
where M and o are positive constants. 
Proof: Let de = wn + Toys kK eV -ells  psenis.be athe 


distinct eigenvalues of A. Since the characteristic nolynomial of A 
is stable, there exists a > 0 such that UE tece- Ul kK = |,. 90m. 
A solution Zz, (t) corresponding to A, is of the form 

7 (3) IS) 
where bp is a constant vector having the entry 1 in one of its rows 
and zeros elsewhere, and where r depends on the multiplicity of ee 
Then since Uy ta < 0, 

|z,(t)|| e%0 as tes 


and thus there exist constants ML such that 


-at 
Zeer es Ber fer >t .. 
If x(t) is any solution to (3.7), it can be expressed as the linear 


combination 


.(t), 


me). = DB. Z. 
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where the summation is over the n basis solutions. Then 
x(t] < 2lsat [l2,(t)| 


< max |6,| = M, are <M aot 


where M = max 18. | =M, . 
It is now easy to show that when the characteristic polynomial 
of A is stable, then (3.7) is an asymptotically stable system. 
Proposition 3.3. 1) ene Characteristic polynomial Ofwe iss cdmilem 
then the system x' = Ax is asymptotically stable. 
Proof: By Lemma 3.2, every solution x(t) of (3.7) satisfies 


||x(t) || < Me 


for appropriate M and a. Since the columns of 
X(t) are solutions, where X(t) is the fundamental matrix of (3.7), 
we have 

HIX(t)]]| <Mne* tot, 
and thus ||X(t)||+0 as te» and the system is asymptotically stable 
by Corollary 1 to Theorem 3.1. 

Thus the stability of a linear homogeneous system with constant 
coefficients depends solely upon the nature of the eigenvalues of A. 
The question now arises as to the stability of a system (3.4) in which 
A(t) is not constant. Proposition 3.1 insures that integrability of 
A(t) is sufficient to guarantee stability in the system. However, this 
is a rather stringent requirement and is of little practical signi- 
ficance in determining the stability of a given system. It is often 
easier to determine stability in a given system by comparing it to a 
known stable system. The simolest case in which this may be done 


successfully is when A(t) is asymntotic to a constant matrix B. 


The following theorem then applies 
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Theorem 3.2. If the characteristic nolynomial of B is stable, 


A(t) is continuous on [0,”), and 
Paes a a Fike 
0 
then the system x' = A(t) x is asymptotically stable. 
Proof: The system x' = A(t) x may be written in the form 
ye= Bx + LACE) - Bl x 
and in this manner may be treated as though it were a non-homogeneous 


system with constant coefficients. The solution x(t30,x,) may then 


be written by Proposition 2.2 and (2.13) 
t 
xc = Xe) Ae eS Kis) A(s) - Bl) x({s) ds, 
0 


where X(t) is the fundamental matrix of the constant system 
x' = B x. By Lemma 3.2., there exist M,oa >0O such that 
1|X(t)| | < Ment and thus 


t 
t 
Jix(t)}] e®” <M [ix ot] + sm [fA(s)-BI | [[x(s)]] eds. 
O 
Now applying Gronwalls‘ Lemma, 
at \ 
HIx(t) || e%* <M [ixgi| exp sm [IA(s)-B]| ds, 
0 


from which 


co 


x(t) | e&* <M [[xof] exo Pm [1A(s)-B| | ds 


aa, 


where K is a convenient constant. Hence ||x(t)| | < Keo and if 


as 


y(t30,y) is any other solution, 
I|x(t) - y(t)]] < 2ke™® +0 as to 
and thus the system x' = A(t) x is asymptotically stable. 

Theorem 3.2. may be modified in the event A(t) does not differ 
much from B, i.e., ||A(t) - B/| is small as follows. 

Theorem 3.3. If the characteristic polynomial of B is stable, 
and ||A(t) - B]| < c on [0,~) for sufficiently small c, then the 
system x' = A(t) x is asymptotically stable. 

Proof: Let x(t30,x,) be a solution of x' = A(t) x. Then xt) 


may be written 
1 
ice » X(t) X, + ; X(t-s)LA(s) - B] x(s) ds, 


where X(t) is the fundamental matrix of x' = Bx. Then proceeding as 


in the proof of Theorem 3.2, we have as before 
at to 
aipciihl-e "<M LXa I | exp Ms |{A(s) - BI] dss 
O 


where M and a are as previously defined in Theorem 3.2. Thus 
L[x(t)[] < M [Ix] | exp (Mc-o) t , 
sO tidl if € 1S Suimielonuly small anieier. C= TE then x(t)-0 as 
teo and since x(t) is an arbitrary solution, the system is 
asymptotically stable. 
If the characteristic polynomial of B is not stable but the solu- 
tions of x' = B xX are bounded, then the following theorem holds. 


Theorem 3.4. If all solutions of x' = B x are bounded on [0,@), 


A(t) is continuous on [0,) and 


ft AGr eel led tac. 
0 
then the system x' = A(t) x is stable. 
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Proof: Let X(t) be the fundamental matrix of the constant 
system x' = B x. Then there exists M > 0 such that ||X(t)|| < ™ 


for teu. IT x(t30,x,) Tseaesoluctiow*oT x’ = gt) x, then 
t 


wit Xe ) xe te X(t-s)[A(s) - B] x(s) ds. 
0 
Thus 
1 

| Pecea) |< Memelx | list Be ats) = B| | |xteenl| wes , 

and applying Gronwall's Lemma, 
t 
[Lx(t)1] <M [Lxgl| exp £m [1ACs) - Bl ds, 


from which 


| |x(t) | | 


| ra 


M |[x [| exp sM |JA(s) - Bl| ds 
0 


Sasecs 


| N 


where K is a suitable constant. The above shows that all solutions are 
bounded and hence the system x' = A(t) x is stable by Theorem 3.1. 
Theorems 3.2, 3.3, and 3.4 insure that if ||A(t) - B|| is well 
behaved, then the solutions of x' = A(t) x will resemble those of the 
constant system x' = B x as tm when A(t) is asymptotic to B. In 
view of this, we might expect these results to generalize to the case 
when B is not constant. This is not always the case though it is 
possible to effect a comparison between two non-constant systems under 
more restrictive hypotheses. Toward this end, we prove the following 


lemma. 


Lemma 3.3. If the system x' = A(t) x is stable on Lt sad 
and 
It 
Lim 7G JA coegAcy edo eer 
t+ i 
0 
then x7 4) is bounded on Lt»), where X(t) is the fundamental matrix 
of the system. 
Proof: Let X(t) be the fundamental matrix of the system. 


Then 


ce = Be 


and since x' = A(t) x is stable, it follows that X(t) is bounded for 
t > t, and this implies that there exists L > 0 such that 
Pragg X(t) | | ae ron t ieee 


Then this relation and (2.8) yield 


xo (t) || = t ; 
exp fs tr Als) ds 
ty 
and since 
t 
eer ae Nees antes 
t+ ty 


it follows that 


exp fs thwA(siads “Suuke> 0’) 


I 
0 


so that 


- | k 
| Neale ila Toueahar t atl 


and thus x7 e+) is bounded. 
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We can now state the following comparison theorem for two non- 
constant systems. 


Theorem 3.5. If the system x' = B(t) x is stable on [t»-), 


ie 
laweinf £5 tr Bis) ds s.- o, 
t->00 iE 
0 
A(t) is continuous on [t 2) and 
s ||A(t) - B(t)|| dt < », 
i 
then the system x' = A(t) x is stable. 


Proof: Let X(t) be the fundamental matrix of the system 


x' = B(t) x. Then the hypotheses imply that there exists M > 0 
such that for t >t. ||X(t)|| <M, and [[X7'(t)|| <M. 
Then if xX(t3t, Xx) is any solution to x' = A(t) x, x(t) may be written 


rc 
x(t) = X(t) xo + s X(t) X7'(s) [A(s) - B(s)]x(s)ds, 


ive 


and utilizing the above we obtain 


ie 
Hx(t)}] <M [xl] + 2 Me [JA(s) = BGs) | []x(s) I Ids. 
tO 
Applying Gronwall's Lemma yields 


t 
HIx(t)|] <M [xpi] exe £m IIA(s) ~ B(5){| ds, 


0 
so that 
[lx(t)|| <M [[xg|| exo sm? ||A(s) - B(s)| [ds 
t 
0 
hin 
where K is a convenient constant. Thus all solutions of x' = A(t) x 


are bounded and hence stable by Theorem 3.1. 


Za 


Theorem 3.5 gives one set of conditions under which a comparison 
may be effected between two non-constant linear homogeneous systems 
and is a well known result. Another set of conditions which have 
been formulated by the author and which are also sufficient to effect 
a comparison are stated below in the following theorem. 

Theorem 3.6. If the system x' = B(t) x is stable on [t.~), 


A(t) is continuous on [t») and 
r qixe! 
it 
0 

where X(t) is the fundamental matrix of x' = B(t) x, then the 


(OME = BCtIeX(t)| | dive» 


system x' = A(t) x is stable on [t»~). 


Proof: Let X(t) be the fundamental matrix of x' = B(t) x. Then 


ie x(t3t) x) is any solution to x' = A(t) x, x(t) may be written 
: 1 
.Gpie= X(B) el KGex «6(s) [A(s)-B(s) ]x(o)iaee 
t 
O 


From this we obtain 


x eyx(t) = x + yp xo! s)LAls) - Ble)]x(s) ds 


and thus 
[xT (t)x(t) | < [xo] + 


oe re 
r | |X” ' (s)[A(s)-B(s)IJX(s)] | | |X" (s) 


Xiewal | ds, 
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so that applying Gronwall‘s Lemma yields 


x7! (e)x(t) | | < [1x1] exp f 1x7! (s)LA(s)-B(s)] 
ms) || ds. 
Therefore, 
DCM edCEYIT x Ixgl | exe [DCMS ) EAs) 8681 
i X(s){| ds 
& K < 
which shows that x7  (t)x(t) is bounded. 
However, 
xe] = UP xCeyx7 (eye (td i] < 1UxCed | 1x7! (t) 
||| 


so that ||x(t)|| < K ||X(t)]] and since x' = B(t) x is stable, X(t) is 
bounded and thus x(t) is bounded. Hence all solutions to x' = A(t) x 
are bounded and thus stable by Theorem 3.1. 

Theorems 3.2 through 3.5 provide conditions under which the 
stability of a given system may be determined by a comnarison with 
a known stable system. Other conaitions under which similar 
theorems can be established have been given by Bellman [1], and 


Levinson [6], as well as numerous other authors. 


ZS, 


TV. NON-HOMOGENEOUS LINEAR SYSTEMS 


The systems considered here are of the form 
(4.1) = Aes Dict 
where A(t) and b(t) are continuous on [t 2). 

Recall that for the homogeneous system, b(t) = 0, boundedness 
and stability of the system are equivalent. This is no longer true 
for the system (4.1) as the following example illustrates. 


Example 4.1. The one-dimensional system 


x' = 2t on [0,-) 
has solutions x(t30,x) ) of the form 
d 2 
x(t) = xX, +t 


which are stable yet unbounded, 

The above example shows that it is possible for a non-homogeneous 
to be stable but unbounded. The reverse behavior, however, is not 
possible as the following theorem shows. 

Theorem 4.1. If the system x' = A(t) x + b(t), A(t) and b(t) 
continuous on [t»2); is bounded on [t»°) then it is stable on 
[t 2°). 

Proof: Let X(t) be the fundamental matrix of the homogeneous 
system x' = A(t) x. If (4.1) is bounded, then in particular, the 
solution x(t3t, »0) given by 


t 
x(t) = X(t) ¢ x's) bls) ds 
it 
O 


is bounded. Thus there exists M)> 0 such that ||x(t)|| <M, t>t. 
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lee y(tst) s¥) is any other solution to (4.1), then 


D 
y(t) = X(t)y + Fs X(t)x7!(s) b(s) ds 
te 
and there exists M., such that 
Hly(t)|] = [|X(t) y, + x(t)|] <M for t > t. 
Then 
x(t) ysll - [Ix(t) || < My 
so that 


eo) yew ata, 
where M = M, + Mo. Therefore if w(tst ow) and v({tst sv) are any 
two solutions of (4.1) and if « > 0, then for t > t, 

FIw(t) = vt)[] = IRC) (wyevy)IL <M LIMg=¥Q I 
so that if ||w-v || < w> then ||w(t)-v(t)|| <e ¢t >t 
which shows that (4.1) is stable. 

Theorem 4.1 admits the following partial converse. 

Theorem 4.2. If the system (4.1) is stable and one solution is 
bounded, then all solutions are bounded. 

Proof: Let X(t) be the fundamental matrix of the homogeneous 
system x' = A(t) x. Once again the first part of the proof consists 
of showing that X(t) is bounded. 

If (4.1) is stable, then in particular, the solution x(t3t 0) 
given by 


t 
x(t) = s X(t)x7!(s)b(s) ds 
tO 
is stable. Therefore, if « > 0, there exists 6 > 0 such that 


fort. > bY, 
aoa O 


3] 


}iw(t)-x(t)|] = [[X(t)w |] < « 


whenever w(t;t ) is a solution with LIwol < 6. Letting 


W 
0° 0 


Woes xe, where e' is the vector whose jth comnonent is e, = §.., we 


1J 
obtain from the above 


where #'(t) is the ith column of X(t). Then since this is valid for 


[|x(t)|| < BE =m, 


Now let z(t3t,»Z)) be a bounded solution and y(tst, sy) any arbitrary 
solution of cae Then since y(t) = z(t) + [y(t) - z(t)], 
ytd] < Eley] + ixct)Ly,-2 111 
<||z(t)|] +™ [ly, = 2oll 

which shows that y(t) is bounded. Hence the system (4.1) is bounded. 

Theorems 4.1 and 4.2 show that in a non-homogeneous linear system, 
boundedness is a stronger property than stability. For this reason, 
it is usually easier to determine stability in a system (4.1) than it 
is boundedness. In fact, the following proposition Shows that the 
stability of (4.1) may be determined by simnly examing its homo- 
geneous part. 

Proposition 4.1. If the system x' = A(t) x is (asymptotically) 
stable on [to~)> then so is the system x' = A(t) x + b(t), where 
bites 1S conti mlolismon [t >). 

Proof: If x(t3t»Xx,) and y(tit,s¥) are two solutions of 
x' = A(t) x + b(t), then the difference [x(t)-y(t)] satisfies the 
homogeneous system x' = A(t) x. Thus since the trivial solution 


z(t) = 0 of the homogeneous system is stable, given c« > 0, there 
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exists 6 > 0 such that if ||(x, - yo) - O|| < 6 then 


0 
WMeatele= Yt) Ope tat, 
which is nrecisely the statement that the non-homogeneous system is 
Stable. The proof is the same if the homogeneous system is, in 
addition, asymptotically stable. 

If A(t) is a constant matrix A, then (4.1) becomes 
(4.2) x' = AxX + b(t) 
and boundedness is easier to determine. If the characteristic 
polynomial of A is stable, then the following theorem holds. 

Theorem 4.3. If the characteristic polynomial of A is stable, 
b(t) is continuous and ||b(t)|| < K, for t > 0, then the system (4.2) 
is bounded. 

Proof: Let X(t) be the fundamental matrix of the constant 
system x' = A x. By Lemma 3.2, there exist M and a such that 


||X(t)| | < Me™ 


le 16 UR x(t3t, x) is a solution of (4.2) 
then from 

t 
x(t) = X(t)x, + s X(t-s)b(s) ds, 

0 


we deduce successively 


a(t-s) 


t 
lix(t)]} <M |[xf| + of Mem K ds 


SM Hy 4. “ (J-e72*) 
<™ I[xt| + 
which shows that x(t) is bounded. Thus the system (4.2) is bounded. 
A similar theorem holds when the solutions of x' = A x are 


bounded, though the forcing term, b(t), must be integrable. 
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Theorem 4.4. If the constant system x' = A x is bounded on 


[0,0), b(t) is continuous on [0,~) and 


Peeiecat) dt erm 
0 


then the system x' = A x + b(t) is bounded. 
Proof: Let X(t) be the fundamental matrix of x' = A x. Then 
there exists M > 0 such that ||X(t)|| <M t>0O. If x(t30,x,) is 


a Soliton tosx! = A. taab( tigithens from 


we deduce successively 


t 
PPA ss Tl PS is: : M ||b(s)|| ds 
< Miasliix | ste | bis i) leds 
- 0 
0 
Sloe 
where K iS a convenient constant. Thus the system x' = A x + b(t) 


1s bounded. 

The boundedness of a non-homogeneous system x' = B(t) x + b(t) 
when B(t) is asymptotic to a constant matrix A may be determined by 
effecting a comparison between it and the constant svstem x' = A x. 
The following theorem gives one instance when this may be done 
successfully. 

Theorem 4.5. If the characteristic polynomial of A is stable, 
b(t) is continuous and bounded on [0,~), B(t) is continuous on [0,-) 
and jf (ttBtt) - A|| dt < », then the system x' = B(t) x + b(t) is 


bounded and asymptotically stable. 
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Proof: Let X(t) be the fundamental matrix of the constant system 


x' = Ax. Then there exist M, a > O such that ||X(t)] | < Met 


font >0.. If x(t30,x,) is a solution to x' = B(t) x + b(t), 
then as before, 
is 
rior = X(t)x, + ys X(t-s)[(B(s)-A)x(s) + b(s)] ds, 
0 
from which 


t 
HIx(t) |] < TIXCe) TT TEx tt + ; | |X(t-s)| | 


(||B(s)-A|| ||x(s)]| + [[b(s)[]) ds. 
Then since ||b(t)|| < N 


t 
[Ix(t) |] e** om [fx] + Me7“S||B(s)-A|| | |x(s) | |Jds 


t 
+ s MNe °° ds 
0 
c 
<M [[xo]] + £Me°]Ix(s) || |1B(s)-Al ds 
+ s MNe*> ds 
0 
and thus 
at -as 
Beso) | ates ee f Me "| |[x(s)|| [|]B(s)-Al] ds, 
where L is a convenient constant. Apnlying Gronwall's Lemma vields 
‘i L 
||x(t)| Jer" < L exp s M |[B(s)-A]| ds, 
0 
so that 


co 


I}x(t){Jet® < Lb exp iM ||B(s)-A|| ds < K< @ 
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where K is another constant. Thus for t > 0, |{x(t) || < ke oe cul 
which shows that x(t) and hence the system x' = B(t) x + b(t) is 
bounded. 

In addition, if y(t30,y,) is any other solution to this system, 

H[x(t) = y(t) i] < [Ix(t)|] + [[y(t) |] < 2ke7e* 

which shows that ||x(t) - y(t)||+0 as t>o and the system 
x' = B(t) x + b(t) is also asymntotically stable. 

If the characteristic polynomial of A is not assumed stable, 
but the system x' = A x is bounded, then a similar theorem holds. 

Theorem 4.6. If the system x‘ = A x is bounded on [0,~), b(t) 


and B(t) are continuous on [0,-@), 


f ||b(s)|] ds < », and 
0 


f||B(s) - A|| ds<o@, 
0 


then the system x' = B(t)x + b(t) is bounded. 
Proof: Let X(t) be the fundamental matrix of the constant 
system x' = A x. Then there exists M > 0 such that ;|X(t)|| < M for 
CeO sl x(t30,x,) is any solution to x' = B(t) x + b(t), then from 
ib 
x(t) = X(t)xX, + f X(t-s)[(B(s)-A)x(s)+b(s)J ds, 


we have successively 


t 
P(e) cM Vix + £m 1B(s)-Al] I]x(s)I| ds 


Ct 
+ s M ||b(s)|| ds 
0 
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He 
eax) > ; M ||B(s)-A]| [|x(s)]] ds 
y , M ||b(s)|| ds 
0 


ie 
K + : M|{B(s)-Al| |[|x(s)]] ds, 


| A 


where K is a suitable constant. Applying Gronwall's Lemma yields 


1 
|x(t)|| <Kexp s M |[B(s)-Al| ds 
0 


< Kexp s M ||B(s)-A[| ds < L < 
0 


where L is another constant. Thus x(t) and hence the system 
x' = B(t) x + b(t) are bounded. 

In general, if A(t) is not constant, Theorems 4.3 through 4.6 
need no longer apply. Thus if a homogeneous system x' = A(t) x 
is bounded (stable) and 

PLPC dt <=, 


it need not follow that the non-homogeneous system x' = A(t) x + b(t) 

1s bounded. The same applies if the homogeneous system is asymptotically 
Stable. A similar problem was encountered in attempting to effect a 
Stability comparison between two non-constant linear homogeneous systems. 
There it was noted that additional hypotheses were required to insure 

the stability of the system x' = B(t) x when B(t) is asymptotic to a 
non-constant matrix A(t) and the system x' = A(t) x is stable. One 


reason for this is that stability, as defined in the sense of Lyapunov, 


3/ 


is quite delicate and often may not be maintained even when small 
perturbations are introduced. This leads to the concepts of uniform 


and restrictive stability, which are discussed in section \V. 
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V. UNIFORM AND RESTRICTIVE STABILITY 


Since Lyapunov stability does not suffice to insure certain 
boundedness and comparison properties, we are led to a consideration 
of stronger forms of stability. One of these is uniform stability, 
defined as follows. 

Definition 5.1)....bet x(t3t sx) be a solution of the system 
x' = £(t,x), where f(t,x) is defined and continuous in the region 
m— 1 (ten: ti <t<, ||x|| < € }. Then x(t) is said to be uniformly 
Seable if it _satisfies (3.21,11, and i113) and im addition (5.1) 
given « > O there exists 6 = &(c,f) > O such that every solution 


y(tst sy) with | |x(t, ) - y(t, )|| < § for some t, > t, satisfies 


eaie= WAT) || ¢< c Fort t,. 

Comparing the definitions, we note that uniform stability is 
Similar to Lyapunov stability. However, there is an important 
distinction between the two. A solution x(t) is stable in the sense 
of Lyapunov if any other solution which is initially close to x(t) 
remains in an e-neighborhood of x(t) for t > ty; at later times, 
me... t > ty} other solutions may come arbitrarily close but need 
not remain there. If x(t) is uniformly stable, once a solution enters 
a prescribed small neighborhood of x(t) for any t, ane jt must remain 
there for all t > t,. Thus, in uniform stability, 6 no longer 
depends on to: 

The above remarks show that a uniformly stable solution is stable. 


If all solutions to a system are uniformly stable, we say that the 


oe, 


system 1s uniformly stable. For linear homogeneous systems of the 
form 
i5.2) Se estar x, 
if one solution to the system is uniformly stable, then the system 
is uniformly stable. Uniform stability for a system (5.2) may be 
characterized by the behavior of the fundamental matrix of the 
system as the following theorem shows. 

Theorem 5.1. The system x' = A(t) x is uniformly stable iff 


there exists M > 0 such that Hx(t)x7 (8) |] <M t> s >t, where 


O 
X(t) is the fundamental matrix of the system. 

Proof: Assume that the system x' = A(t) x is uniformly stable, 
so that the trivial solution z(t) = 0 is uniformly stable. Then 
if « > 0 there exists 6 > O such that if | |x(t,) | | < § for some 
if 


Seta. Chen. | Cty eee: TOR te > ty. Let y(t3t sy) be a solution 


| Oo: 
of x' = A(t) x. Then y(t) may be written 
y(t) = x(t) x7 '(t,) ylt,). 


By choosing y(t) to be that solution satisfying y(t.) = se, where 
e is as in the proof of Theorem 3.1, we have | ly(t,) | | = >< 8 


so that 


Sin ( DIAS bistede lilo: me aemt t 


or equivalently, 


Ze 
ace >) 





i (exe Mic 


where alt stew) is the ith column of the matrix cae Repeating 
the above procedure for i = 1,...,n yields 


-| en 
PPG) Ce eS ee 
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Finally, since this procedure is valid for any ty ida and the 
bound is indevendent of ty» we obtain 
-| 2n 
POM Ge <a set potpeeet 
Conversely, if ||X(t) X7'(t,)|| <M for t> t 





9° 
wi to: and 
x(t3t, »X,) and y(tst,s¥,) are any two solutions of x' = A(t) x, 


then for t > ty ei to 


| |x(t)-y(t)| | 


|X (t)x7! (t, Ex(t, )-y(t, I] | 
M | |x(ty) - y(t, 1 


so that if ec > 0, ||x(t,) - y(t))|| <q t) >t, implies 


| A 


Leroy yc || <2 t > ty which shows that the system is uniformly 
stable. 

From Theorem 5.1, we obtain the following simple proposition. 

Proposition 5.1. If the system x' = A x, where A is constant, 
is stable on [0,-), then it is uniformly stable on [0,~<). 

Proof: Let X(t) be the fundamental matrix of x' = A x. Then 
there exists M > O such”that ||/X(t)|| < M for t > 0. Then by (2.13), 

HixX(t-s)}| = []x(t) x7 '(s)|| <M for t>s>0 
and the system is thus uniformly stable. 

Recall that in sections III and IV, several boundedness and 
comparison properties which hold when A(t) = A is constant, need not 
apply when A(t) is not constant. However, Proposition 5.1 shows 
that in the case A(t) is constant, the systems considered are actually 
uniformly stable, rather than just stable. This leads to the question 
of whether or not uniform stability of the non-constant svstem 
x' = A(t) x is sufficient to insure these boundedness and comnarison 
properties. The following theorems provide an affirmative answer 


to this question. 


4] 


Theorem 5.2. If the system x' = A(t) x is uniformly stable on 


[t 2). B(t) is continuous on [t.»~) and pf | [B(t)-A(t) | | dt 
O 


then the system x' = B(t) x is uniformly stable. 

Proof: Let X(t) be the fundamental matrix of the system 
x' = A(t) x. Then ||X(t) XH os)] aM tor tes > te leg eee O 
and fix s >t. If x(t3t>x)) and y(tst, sy) are two solutions of 
the system x' = B(t)x, then since B(t) = A(t) + [B(t) - A(t)], 


we nave TOY ts. 


x(t)-y(t) = x(t)x7'(s)Ex(s)-y(s)] 
it 
+ op X(t)X7 (2) [B(r)-ACt) Ex (2 )-y(t)] de 
S 
and thus 


Hix(t)-y(t) || < [[xCt)x7'(s)] | | [x(s)-y(s) || 


t 
ts | x(t) Med] | [BCe)-ACz) | | | x(2)- 
S 


y(t)| [dr 
< M [|x(s)-y(s)|| 


1 
= x M ||B(r)-A(z)|| | |xGz)-y(r)]| dt. 
O 


Applying Gronwall's Lemma yields 
t 
| |x(t)-y(t)| | < M ||x(s)-y(s)|| exp i M ||B(r)-A(x) | | 
. di 


from which 


Hx(t)-y(t)| |< Me] |x(s)-y(s) |] exp ry /B(x)- 


NCOP de. 
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Therefore, 
| |x(t)-y(t)]| < L | |x(s)-y(s) || 
where L is a convenient constant. The above relation is valid 


E 


for t>s>t,, so that if ||x(s)-y(s)|| <P $s >t,» then 


O 
| |x(t)-y(t)|| < © for t > s which shows that the system x' = B(t) x 
is uniformly stable. 

Corollary 1. Under the same hypotheses as Theorem 5.2, the 
system x' = B(t) x + b(t), where b(t) is continuous on [t,») is 
uniformly stable on [t, 2°). 

Proof: The proof follows from the fact that under the hynotheses, 
the system x' = B(t) x is uniformly stable and that if x(t) and 
y(t) are two solutions to the non-homogeneous system x‘ = B(t) x + b(t), 
then the difference x(t) - y(t) is a solution to the non-homogeneous 
system x' = B(t) x. 

When the homogeneous system is uniformly stable, then 
boundedness of the non-homogeneous system may be determined as 
for the constant case considered earlier in Section IV. 

Theorem 5.3. If the system x' = A(t) x is uniformly stable 
on [t om) b(t) is continuous on [t, 2°) and ft Hbts) TI ds < , 
then the system x' = A(t) x + b(t) is bounded on [t,»°). 

Proof: Let x(t3t Xx) be a solution to x' = A(t) x + b(t). 
Then 

ie 


CEM 4 NE) Kes )b(s Mus; 
tO 
where X(t) is the fundamental matrix of the homogeneous svstem 


arty x. Since ||X(t) x7 sy] Samo toms > L, 
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t 
CENT [xl] + fH | Ib(S) I ds 


0 
and thus 
JIx(t) || <M [ix |] + 7M |] b(s)]] ds < b< @ 
be 
where L is a convenient constant. Thus the system x' = A(t) x + b(t) 
is bounded. 


It is also possible to determine boundedness of a non-homoqeneous 
system through comparison with a uniformly stable system. 
Proposition 5.2. If the svstem x' = A(t) x is uniformly stable 


- on [t.~)> B(t) and b(t) are continuous on [t,»~); Ay Gan dt < «, 
% 0 
and s||B(t) - A(t)|| dt < », then the system x' = B(t) x + b(t) 
ty 
is bounded on Wee) 


Proof: By Theorem 5.2, the hynotheses imply that the homo- 
geneous system x' = B(t) x is uniformly stable on [t,»»)- The 
result then follows from Theorem 5.3. 

Thus certain boundedness and comparison properties established 
for constant systems may be generalized to non-constant systems 
which are uniformly stable. 

We now introduce the concept of restrictive stability. 

Definition 5.2. The linear homogeneous system x' = A(t) x 
is said to be restrictively Stable if it is stable and its adjoint 


i x is stable. 


system x = flee) 
We note that restrictive stability applies to a system, rather 


than a particular solution and is defined for linear homogeneous 


4a 


systems only. The following theorem characterizes restrictive 
stability in terms of the behavior of the fundamental matrix of the 
system. 

Theorem 5.4. The system x' = A(t) x is restrictively stable 
on ERS) iff there exists M > 0 such that ||X(t)|| < M and 
ley ils M foget >t» where X(t) is the fundamental matrix 
of the svstem x' = A(t)x. 

Proof: The result follows by Theorem 3.1, noting that rx eeyq! 
is the fundamental matrix of the adjoint system x' = -A(t)'x. 

From Theorem 5.4, we see that a restrictively stable system 
x' = A(t) x is uniformly stable since 

x(x (sd < TICE TDC) <M? sate t, 

but not necessarily conversely. Also, since -f-a(t)!y! =e), 
it follows that if x' = A(t) x is restrictively stable, then its 
adjoint system is restrictively stable. The following propositions 
give conditions which suffice for a system to be restrictively 
Stable. 


co 


Provosition 5.3. If s ||A(t)]]| dt < ©, and A(t) is continuous 
’ 
on m=)! then the system x' = A(t) x is restrictively stable on 


[t 2). 

Proof: Proposition 3.2 shows that under the hypotheses, the 
system x' = A(t) x is stable. Since ||A(t)|| = || -A(t)! ||, it 
follows that the adjoint system is stable. 

Proposition 5.4. If the system x' = A(t) x is stable on [t s) 


and A(t) is continuous on [t2)> and lim inf ice trWAts) dsv> -@, 
T+ t 
0 


then it is restrictively stable. 
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Proof: Let X(t) be tne fundamental matrix of x' = A(t) x. Then 
X(t) is bounded for t > t, since the system is stable. By Lemma 3.3, 


Asie ita f tr A(s) ds > -@ implies that x7! (+) is bounded for 
to 1E 


O 
iC > t. Hence the system x' = A(t) x is restrictively stable by 
Theorem 5.4. 
Proposition 5.5. If the adjoint system x' = -A(t) "x is stable 


on [t»~); A(t) is continuous on Lt»); and lim sup r* tr A(s) ds < , 
t+ t 


O 
then the system x' = A(t) x is restrictively stable on [t»@). 
Proof: Let X(t) be the fundamental matrix of x' = A(t) x. Then 
seers adjoint system is stable, x7! (¢) is bounded for t > t. 
From the relation X(t) = det X(t) adj xTeey, we obtain by (2.8) 


t 
X(t) = exp s tr A(s) ds adj X7!(t) 
t 


O 
and since aes is bounded, adj x Te) is bounded and thus 
t 
}|X(t)|] <Mexp f£ tr A(s) ds 
0 


c 
<Mlim sup s tr A(s) ds <L<o 
t+ ie - 


O 

where M and L are suitable constants. Hence X(t) is bounded and 

the system x' = A(t) x is restrictively stable by Theorem 5.4. 
An interesting result concerning the well known second order 

Hill's Equation 

(5.3) u'' + a(t).u = Om on [Oge)), 


where a(t) is real valued and continuous on [0,-) andeu(t)eais 


a scalar valued function, may be obtained through the use of 
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stability theory. The vector analogue of (5.3) is the two-dimensional 
system 

(5.4) y' = Zz 

a(t) y 

which is in the form (5.2) with x =(¥%) and A(t) = 


7! 
] 


~a(t) 0 | 
We note that, since tr A(t) = 0, if the system (5.4) is stable, 
then it will be restrictively stable, and thus uniformly stable. 


It is now possible to give the following result: If fsfa(t)| dt < , 
0 


then there exist unbounded solutions of (5.3). This may be 

proved using stability theory as follows. Assume that all solutions 
to (5.3) are bounded. Then (5.4) is a stable system and hence 
uniformly stable by the above remarks. Theorem 5.2 then asserts 


Ena wy 


then the system x' = B(t) x is uniformly stable and hence bounded 


co 


since fs ||B(t) - A(t)|| dt = s Ja(t)| dt < ». But this is a 
0 0 


contradiction since 


0 (3 


is a solution to x' = B(t)x. Thus all solutions to (5.3) are bounded. 
Restrictive stability may also be characterized in terms of 
the reducibility of the system. Toward this end, we have the 


following definition. 
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Definition 5.3. Tne system x' = A(t) x is said to be reducible 
if there exists a matrix L(t) whose elements are differentiable 


functions which is bounded together with aes on [t,»~) such 


that 
(5,5) uve) att) L(t) - UT (t) L'(t) = B 
where B is a constant matrix. If B = O, the zero matrix, then 


we say that the system is reducible to zero. 

The relation (5.5) is often exnressed by saying that A(t) is 
kinematically similar to the constant matrix B. If both L(t) and 
A(t) are constant, then (5.5) becomes 


ER i 


which is the ordinary similarity relation. 

If a linear homogeneous system is reducible, then the trans- 
formation x = L(t) y reduces the system x' = A(t) x to the constant 
system y' = B y, while if x' = A(t) x is reducible to zero, then 
this transformation yields the system vy’ = 0. 

The following theorem characterizes restrictive stability in 
terms of the reducibility of the system. 

Theorem 5.5. The system x' = A(t) x is restrictivelv stable on 
[t,»~) iff it is reducible to zero on [t»@). 

Proof: If x' = A(t) x is restrictively stable, then there 
exists M > 0 such that ||X(t)|| < M and iix7 (ty |] <M fort>t, 
where X(t) is the fundamental matrix of the svstem. Then since 


X(t) satisfies 0 t)e= 9A TT tefoMows *that 


xe) alt) x(t) - x7 ey xt) = 0 


and that the system x' = A(t) x is reducible to zero by taking 
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Conversely, assume that (5.5) is satisfied for B = 0. Then 

from (5.5), L(t) satisfies 
L' = AK®) 1 

and since L(t) is non-singular, it follows that L(t) is a basis 
matrix for the system x' = A(t) x. Hence X(t) = L(t) W, where X(t) 
is the fundamental matrix and W is a constant non-singular matrix. 
Thus ||X(t){| < [LCE | WI} and pxe'eyyy < yum tceyy] (wey which 
shows that both X(t) and xl) are bounded. Thus the system 
x' = A(t) x is restrictively stable by Theorem 5.4. 

Uniform stability is also related to reducibility as the 
following theorem illustrates. 

Theorem 5.6. If the system x' = A(t) x is stable and reducible 
on [0,-), then it is uniformly stable on [0,-). 

Proof: Let X(t) be the fundamental matrix of x' = A(t) x and 
let L(t) be as in Definition 5.3. Then the transformation 
y(t) = ues x(t) reduces x' = A(t) x to the svstem y' = B y, 
where B is constant. Also, since eles X(t) is non-singular and 
satisfies 

(eee sr x), 

it follows that L(t) X(t) is a basis matrix for y' = B y. Then 
Y(t) = ees) X(t) W, where Y(t) is the fundamental matrix of 
y' = B y and W is a constant non-singular matrix. Since x' = A(t) x 
is stable, X(t) is bounded and by definition, both L(t) and un! ¢t) 
are bounded. Thus there exists M > 0 such that ||X(t)|| <M, 
[|L(t)| | <M, and |[L7'(t){] <M, so that] |y(t)]| < Mo [|W] |, which 


shows that Y(t) is bounded and thus that the system y' = B y is stable. 
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Then by Proposition 5.1, y' = B y is uniformly stable and thus there 


exists N > 0 such that ||Y(t) v7 (s) J < Ne fort > si> 0. Then 


Since 


X(t) = L(t) v(t) wT! 


we have for t > s > 0 


Pixceyx7 (sd ty < [LCen | Hyeyy7 sq] eel ce) | 


< MON 


| A 


and thus x' = A(t) x is uniformly stable by Theorem 5.1. 

The above remarks relate the connection between uniform and 
restrictive stability in a system and the reducibility of that 
system. Also, we have seen that uniformly stable systems and 
restrictively stable systems satisfy certain boundedness and 
comparison properties that need not necessarily hold for stable 
systems. In Section VI, we shall show that a similar behavior 


is exhibited when the systems considered are non-linear. 
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VI. NONLINEAR SYSTEMS 


The systems considered here are non-linear mixed systems of 
the form 
(6.0) al Ni el ef (ut cel 
where A(t) is continuous on [0,~) and f(t,x) satisfies certain 
appropriate conditions. Such systems are extremely important 
Since the general first order normal system 
(6.2) Mel Gilat Xs) 
can often be put in the form (6.1) by taking f(t,x) = a(t,x) - A(t) 
where A(t) depends on the form of g(t,x). Conditions on f(t,x) 
are thus sought which insure stability in the non-linear system 
(6.1) when it is known that the linear portion 


(6.3) Kus 


A(t) x 
is stable. When A(t) = A is constant, then (6.1) becomes 
(6.4) <! 


Ax + f(t,x) 
and stability is easier to determine. 
Consider the following conditions on f(t,x). 
(6.51) f(t,x) is continuous in the region 
B= i(t,x): 0-< t < =, |{x}} < 2} 
(6.577) L}f(t.x)|] = o(||x|]) uniformly in t as 
||x||>+0, i.e., there exists a function u(r) > 0 
such that u(r)>0 as r->0 and 


H]FCtsx) TT < [ixt] u€{ xf]. 


5 | 


Condition (6.51) assures existence of solutions for small 
||x|| but not uniqueness. Note also that (6.5ii) implies that 
Z(t) = Os a solution to (Gee). 

When f(t,x) satisfies (6.51,i11), we might heuristically 
reason that since f(t,x) is small when ||x|| is near zero, the 
system (6.4) might be expected to behave in a manner similar 
to x' = A x. Such reasoning is correct as the following theorem 
illustrates. 

Theorem 6.1. If the characteristic polynomial of A is stable 
and f(t,x) satisfies (6.51, ii) then the solution z(t) = 0 of 
the system x' = A x + f(t,x) is asymptotically stable. 

Proof: We first show that solutions x(t30,x,) exist and are 
defined on [0,~) when Hx is small. Let X(t) be the fundamental 
matrix of the constant system x' = Ax. Then X(t) satisfies for t > 0 
[|X(t)|] < ™ et for suitable M and a. The solution x(t30,x.) of 
x' = A x + FCCSX) satisfies (where it exists) 

ie 
cy) = ee x, + f X(t-s) f(s,x(s)) ds, 
so that 


ie 
[[x(t)|fe% <M [|x fl] + £ Me*S||F(s,x(s))|| ds. 
0 


The above inequality is valid for t [0,t,) where ||x(t)|| < efor 
t [0,t,), assuming LIXoI < &. Then from (6.511), given « > 0, there 
exists 6 > 0 SUCK@Uhat Tome 07 
H1F(tsx) [1 < el [x] 


whenever ||x|| < 6. 


oc 


If we take Hx ot | < 6, then by continuity of x(t), there exists 


Lao Oesuch ‘that™| ject) | |. Fommue<"t < ty. Therefore for 


Z 
0 < t < T, where T = min (t sto); we have 


, 
[Ix(t)|| eF sm [[xgi] + Me se] [x(s) {| as. 
Applying Gronwall's Lemma to the above yields 
meat 
||x(t)| Je < MI [x] | exp sf Me ds 
0 


ang tamis for 0 < t 8h 
[|x(t)]| <M Hx exo(Me-a)t. 


Choosing e < W yields 


HIx(t)]] <M [ix ff] O<t<T 


so that if ||x || < ay > 


lae( tdilale < >< $ Occee < I. 


Since f(t,x) is defined and continuous for t > 0 and ||x|| < &, 
the above relation imolies that tne solution x(t30,x,) can be 
extended interval by interval preserving the above bound. Thus, 


any solution x(t30,x,) with HIxoH < mt exists on [0,~) and satisfies 


Pest tell = 5 for t > 0. Since 6 can be taken arbitrarily small, 
this implies that the solution z(t) = 0 of x' =A x + f(t,x) 
is stable. Furthermore, 
Pee ie Hx oI | exp (Me-a)t 
where © <q implies that ||x(t)||+0 as te» and thus the trivial 


solution js asymptotically stable. 


oS 


Coddington and Levinson [3] show that the hypotheses of Theorem 6.1 
may be weakened somewhat. Instead of requiring that f(t,x) satisfy 
(6.511) it is sufficient to assume instead that for some K > 0, 
HF(tax)]] < K [|x|] t > 0 

for all small ||x||, and that qiven e > O there exist 6,T such that 
#(tsx)I1 < © TIxl| 

for Wix|))< éuendetee. 1. 

Theorem 6.1 may then be applied to a system of the form 
(6.6) x' = Ax + Bit) x + g(t.x) 
where B(t) +0 as te» and ||g(t,x)|| = o(]||x]]). 

Other conditions on f(t,x) which suffice for theorems similar 
to Theorem 6.1 may be found in [1] and [6]. 

When A(t) is not constant in (6.1), then Theorem 6.1 need no 
longer necessarily hold. However, if the linear portion xt = A(t) x 
of (6.1) is uniformly stable, then we have the following theorem. 

Theorem 6.2. If the linear system x' = A(t) x is uniformly stable 


and f(t,x) satisfies 


(6.77) f(t,x) defined and continuous for t > 0, 
||x!| <€ and 
(Geaiip) | f(t,x)|]| < a(t) ||x]], where a(t) > 0 and 


co 


f a(t) dt < 2, 
0 


then the solution z(t) = 0 of x' = A(t) x + f(t,x) ts untformly 
stable on [0,~). 

Proof: Let X(t) be the fundamental matrix of the linear system 
x" = A(t) x. Then there exists M> 0 such that ||x(t) x7!(s)|] <™ 
for t > s > 0. As in the proof of Theorem 6.1, it can be shown that 
any solution x(t30,x,) of x' = A(t) x + f(t,x) exists and is defined 
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on [0,-). Then assuming this, the solution x(t30,x,) can be written 


t 
x(t) = X(t)x7'(s)x(s) + te X(t)X7 (2) #(2,x(r)) dr, 


where s > 0 is fixed, so that 


t 
I]x(t)]] < MI ]x(s)]] + ~ a(x) ||x(+)| | dr. 


Applying Gronwall's Lemma yields 
t 
Iix(t)] |< M[[x(s)]] exp os Ma(r) det 
0 


and thus 


||x(t)]] < M||x(s)|] exp a a(t) dr < K{|x(s)1| 
where K is a convenient constant. The above relation is valid for any 
s > 0 so that given ce > 0, ||x(s)|| < ra for some s > 0 implies 
||x(t)|| < © for all t > s so that the trivial solution z(t) = 0 of 
x' = A(t) x + f(t,x) ts uniformly stable. 

Corollary 1. Under the same hypothesis as Theorem 6.2, if the 
linear system x' = A(t) x is in addition asymptotically stable, then 
so is the trivial solution z(t) = 0 of x' = A(t) x + f(t,x). 

Proof: The hypotheses imply that z(t) = 0 is a uniformly stable 
solution of x' = A(t) x + f(t,x). If the linear system is in addition 
asymptotically stable, then ||X(t)|| + 0 as t+» where X(t) is the 
fundamental matrix of the linear system. Therefore, if x(t30,x,) is any 
solution to x' = A(t) x + f(t,x), then aiven c>0 there exists T such 


that ||X(t){| <= for t>T. Then from 
0 


t 
x(t) = X(t)x, + J x(t)x7 (5) #(s,x(s))ds, 


we have 


t 
Hix(t)]] < [X(t] [ix ot] + - Ma(s) ||x(s)[] ds, 
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where | |X(t) x Ts) <M for t > s > 0. Then for t > T we obtain 


| L 
pleas) | oeeeh 4 a(s) ||x(s)|| ds 
and applying Gronwall's Lemma, 


t 
||x(t)| | < © exp RS M a(s) ds, 


so that eGo) iim xD r Ma(s) ds < eK 

where K is a suitable constant. Since e« >0 1s arbitrary and K is 
independent of « and T, we conclude that ||x(t)||+0 as te» and 
that the trivial solution z(t) = 0 of x' = A(t) x + f(t,x) jis 
uniformly and asymptotically stable. 

From the proof of Theorem 6.1, we note that when the linear 
system x' = A(t) x is uniformly stable and f(t,x) satisfies 
(6.7i1,i1), then all solutions of (6.1) remain bounded. However, 
in the event the linear system is restrictively stable, then the 
solutions of (6.1) are actually bounded by those of the linear 
system as the following theorem illustrates. 

Theorem 6.3. If the system x' = A(t) x is restrictively stable 
and f(t,x) satisfies (6.7i,ii), then every solution x(t30,x,) 

Gfx = NUL) Xt i t.x) Samistiesmauny = OC] | X(t) 1 |). ame... 
PRG) Peers Lx ol where K > 0 is constant. 

Proof: Let X(t) be the fundamental matrix of the linear system 
x' = A(t) x. Then there exists M > 0 such that ||X(t)|| < M and 
Themcall <M fort>0O. Let x(t30,x,) be a solution of 


x' = A(t) x + f(t.x). Then from 


t 
x(t) = X(t)x, 4 J x(t) x7 '(s) #(s.x(s}) Gs, 
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we obtain 
t 
x(t) 
Lett < Lixgl + MF als) [Ix(s) || ds 


et 
2 x(s) 
Sue + M ry a(s) Ye ds. 


Applying Gronwall's Lemma yields 


L 
8 < [ix [| exp s Me a(S) als 
: 0 
M d 
< Ix it exp g a(s) ds <Klix |, 
where K is a convenient constant. Thus | x(t)) || SK ix IIycty | 


which completes the proof. 

Theorem 6.3 shows that when the linear system is restrictively 
stable, the contribuition of the non-linear term f(t,x) to the solutions 
of (6.1) is small. Golomb [4] has given conditions under which bounded- 
ness and stability may be determined in systems of the form in which 


i teexe) Sits fines 


(6.87) Tees iseconbinuous’in X a.e. fort > 0 


and is intearable in t for fixed x and 


(6.8ii) lF(t.x){] < a(t) u€||x]]) for t > 0, where 
a(t) ts positive and intearable and u(r) 
defined for r> 0 is a positive non-decreasina 
continuous function which is submultiplica- 
DIVE. VGar 


u(r s) = uly) Uls), 


Sy 


without requiring uniform stability in the linear system, essentially 
qeneralizing Theorems 6.2 and 6.3. Other evaluations of solutions to 


(6.1) may be found in the same paper. 
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VIT. WUNGEUSTONS 


We have considered stability in the sense of Lyapunov for 
systems of ordinary differential equations with particular emnhasis 
on linear systems. The theory for these systems is quite advanced 
today and in Sections III, IV, and V we have considered several of 
the more important results. Though by no means exhaustive, much of 
the material in these sections comes from varied sources and it is 
felt that the presentation there offers a more complete treatment 
than is usually found elsewhere. In addition, a new set of sufficient 
conditions for a stability comparison theorem have been formulated 
which to the best of the author's knowledge have not avneared else- 
where. Though the treatment in Section VI is somewhat brief, it is 
intended for the most part to complement the rest of the paper by 
illustrating some of the more useful properties exhibited by uniformly 
and restrictively stable systems. The theory of stability for non- 
linear non-autonomous systems is still quite incomplete and is today 
an area of extensive active research. 

Some other aspects of stability theory not considered here are 
the Poincare-Bendixson theory for autonomous systems, the second 
method of Lyapunov, and various analvtical approaches such as the 
method of Lindstedt, Krylov and Bogolyubov, and of Van Der Pol. 

As mentioned earlier, the approach to the problem of stability of 
a system of differential equations varies with both the form and 


order of the system. This paper has considered but one of these 


Se 


approaches, one which is particularly useful in the treatment of 
linear systems. For a treatment of some of the other areas in the 
theory of stability of ordinary differential equations, we refer the 
reader to several of the excellent texts on stability, including 


[24 TS rand" 7" 
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